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Abstract 

, 

We consider the intensity pattern, generated by a monochromatic source, 
in a disordered cavity coupled to the environment. For weak coupling, and 
when the source frequency is tuned to a resonance, the intensity distri- 
bution P(I) is close to Porter-Thomas distribution. When the coupling 
increases, P(I) gradually crosses over to the Rayleigh distribution. The 
joint probability distribution for intensities at two different points is also 
^ ' discussed. 



PACS numbers: 42.20-y, 71.55Jv 



A wave propagating in a random medium produces a complicated, highly ir- 
regular intensity pattern. This pattern is described in statistical terms. One con- 
siders an ensemble of different realizations of the random medium and inquires, 
e.g., about the probability distribution P(T) of the wave intensity I at some point 
r. In many cases P(I) is accurately described by the Rayleigh distribution 

P B (/) = ^yexp(-//(/» , (1) 
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where (I) is the average intensity. A simple derivation of Eq. (1) is based on 
an assumption that the field at point r can be viewed as a sum of many random 



lating corrections to the Rayleigh distribution || ||, is based on the perturbative 
diagrammatic technique. Here one starts with the wave equation 



supplemented by appropriate sources. In this equation ipuij) describes a field (for 
instance, pressure field in an acoustic wave or a component of the electric field 
in an electromagnetic wave), excited by a monochromatic source of frequency uo. 
The random function /z(r) describes the fluctuating part of the refraction index 
and k a = u/c , c being the the speed of propagation in the average medium. In 
the diagrammatic approach one computes moments of the intensity I = |^(r)| 2 
and reconstructs the distribution P(I). 

Let us stress that Eq. (1) applies to the case of a monochromatic wave propa- 
gating in an open system. A different type of problem arises if one considers the 
wave equation (2) in a closed geometry without sources. In this case one inquires 
about the statistical properties of a single eigenstate ip a (r) , e.g., about the dis- 
tribution P{u) of the quantity u = \ip a (r)\ 2 . Extensive studies of chaotic |J and 
disordered cavities have demonstrated that the main part of the distribution 
is described by the Porter-Thomas statistics 



contributions |T], 0. A more systematic derivation, which paves the way for calcu 



V 2 + ^(l + //(r)) ^ w (r) = 



(2) 
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where V is the volume of the cavity and (u) = 1 /V . 

The above discussion suggests the following idea. Let us assume that the cav- 
ity is weakly coupled to the environment, for instance, via a small opening or due 
to some small absorption in the bulk. Then, by placing a monochromatic source 
inside the cavity and tuning its frequency to to a resonance, one will generate in 
the cavity an intensity pattern I(r) which closely follows the "profile" |7/> a (r)| 2 of 
the eigenstate a with frequency u a ~ u>. Thus, the intensity distribution will be 
given by the Porter-Thomas statistics. On the other hand, for a sufficiently strong 
coupling (an open system) the distribution should obey the Rayleigh statistics, 
Eq. (1). The main purpose of this letter is to investigate the crossover between 
weak and strong coupling and to propose a generalized distribution P(I) which 
interpolates between the two regimes. 

Let us first emphasize the difference between an open and a closed system, 
using the simple picture of addition of many random waves 0. In an open system 
the local field ip(r, t) can be viewed as a sum of great number of traveling waves, 
arriving to a point r from various scattering processes: 

N 

?P(r,t)=N- 1 / 2 J2cos(9 n + k n -r-ut) , (4) 

71=1 

where the phases 9 n are completely random and all the amplitudes have been 
taken to be equal (one could assume random independent amplitudes, without 
any change in the results). The wave vectors k n are uniformly distributed on a 
(^-dimensional sphere (d = 2, 3) of radius k Q . The instantaneous local intensity 
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is defined as vp 2 (r,t). The measured quantity, /, is the intensity averaged over 
time, i.e. over one period T = 2tt/u : 

I r T i N 

I=- dt V> 2 (r, t) = — ^s[9 n - m + (k n - k m ) • r] . (5) 
1 Jo 2M n m 

Note that the same expression for I is obtained if one assumes a complex, time 
independent field 

N 

^(r) = (2iV)" 1 / 2 ]T exp[i(0 n + k n • r)] (6) 

n=l 

and defines the intensity as I(r) = ^(r)) 2 . It follows now from the central limit 
theorem that both ^R,ip and Qip are independent Gaussian variables, with zero 
mean and equal variances, which leads to Eq. (1) for the intensity distribution. 
In a closed system the field is viewed as a sum of many standing waves: 

N 

ip(r, t) = N~ l/2 cos(9 n + k n • r) cos cut (7) 

n=l 

As far as the time-averaged intensity is concerned, one can ignore the factor 
cos out and define a stationary field 

N 

^(r) = (2iV)- 1 / 2 ^cos(^ + k n -r) . (8) 

n=l 

The central limit theorem now tells us that i/) u (r) is a Gaussian variable with 
zero mean, and the Porter-Thomas statistics 

/ j xi/2 

Pp - t{I) = [MI)!) ex P(- J /2( J )) ( 9 ) 
for the intensity / = ^> 2 (r) follows immediately. 



After clarifying the difference between sums of random traveling waves (open 
systems) and standing waves (closed systems), we move to the general case of a 
disordered cavity coupled to the external world. We write the stationary field as 



N 



V(r) = [2N(1 + 2e 2 )}- 1 ' 2 £ {cos(#„ + k n • r) + eexp[i(9' n + k„ • r)]} (10) 



n=l 



where the parameter e describes the strength of the coupling. For small e, and 
when the source frequency is tuned to a resonance, the field consists of a large- 
amplitude standing wave (an eigenstate) with a small traveling wave "riding" on 
top of it. The intensity distribution P(I) is close to the expression in Eq. (9). 
Large e corresponds to an open system, where the field is mostly a traveling wave, 
and P(I) is close to the Rayleigh distribution, Eq. (1). 

All the phases, 6 n and 9' n , in Eq. (10) are independent and uniformly dis- 
tributed between and 2tt. It is then clear that both dtif) and ^sip are independent 
Gaussian variables with variances {($ti(>) 2 ) / {(Ssijj) 2 ) = (1 + e 2 )/e 2 . This leads to 
the following distribution for the intensity / = (5fa/>) 2 + (Qip) 2 : 



V ; 2(/)ev / TT^ 



l + 2e 



2\2 



/ 1 + 2e 2 
4(/)e 2 1 + e 2 
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4(/)e 2 1 + e 2 

where J G [x] is the modified Bessel function. The n'th moment of this distribution 
is then given by (I n ) = (I) n n\ 2 Fi [-§ , 1 ^,i,(i+2e 2 )- 2 ], where 2 Fi is the Gaussian 
hypergeometric function. 

Since the parameter e in Eq. (10) is attached to the propagating part of the 
field, it can be related to the (averaged over time) current density0: 

zc 

J(r) = — [^(r)V^(r)-c.c] . (12) 



Substituting tp from Eq. (10) and averaging over phases, one finds that (J(r)) 
vanishes[|] and (J 2 (r)) = 2c 2 (/) 2 e 2 (l + e 2 )/(l + 2e 2 ) 2 . Therefore, instead of using 
the somewhat vague notion of the "coupling strength e" for parametrization of 
the distribution P(I), one can use the dimensionless ratio 5 = (J 2 )/c 2 (I) 2 . 

We, thus, propose a one-parameter distribution Ps(I) for the intensity (we 
normalize I to its average value, i.e. choose (I) = 1): 



The parameter 5 can assume values from (closed system, no current) to 1/2 
(open system, maximal current density). When this parameter changes from 
to 1/2, the intensity distribution changes from Porter-Thomas to Rayleigh. 

Let us mention that a somewhat different crossover phenomenon has been- 
considered in Ref.||. These authors discussed the statistics of |-?/> Q (r)| 2 = u for an 
electron's eigenstate in a quantum dot, in the presence of an arbitrary magnetic 
field. For zero field the distribution P(u) is given by Eq. (3), whereas for suffi- 
ciently strong field it crosses over to a Rayleigh distribution P{u) = V exp(-Vu) . 
In this crossover problem, as opposed to the one considered in the present paper, 
the system always remains closed. 

In a similar way one can consider the distribution for the local current density 
J(r) or the joint probability distribution P(I,J). We will not discuss here these 
objects but limit the discussion to the joint probability distribution P(Ii,I 2 ), 
where = /(r^) (i = 1,2). We start with a wave propagating in an open 
system. The field ip(r) is then given by Eq. (6). It follows from that equation 




(13) 
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that, in the Large N limit, the joint probability distribution for ^(rj.) = ipi an d 
V»(r 2 ) = 4>2 is: 

= ^k^-^UK- 1 )^} (14) 

where Kij = (ipiip*) is the 2x2 covariance matrix with K\\ = K22 = l,i^i2 = 
K21 = f{p), and p = |ri — r2|. The explicit form of the field-field corre- 
lation function f(p) will be given below. Transforming to polar coordinates, 
4>i = V^i ex P(^i)) an d integrating out the phases, one obtains: 



Eq. (14) is the standard assumption in the theory of optical and acoustical speck- 
les and the resulting distribution P(J 1 ,/ 2 ) is well known in optics and acoustics 
of disordered media |], ||. 

After averaging the product ^(ri)^*^) over the random phases 9 n , one finds 
f(p) — Sn ex P(^k n ■ p) where the sum is taken over iV points on a unit sphere. 
In the N —>■ 00 limit, replacing the sum by an integral, one finds f(p) = J (k p) 
in 2 dimensions and f(p) = (A; p) _1 sin(/c p) in 3 dimensions. (A more rigorous 
calculation shows that, for an open geometry, f(p) decays exponentially for p 
larger than the mean free path £). 

Eqs. (14), (15) describe the statistics of radiation in an open system. In con- 
trast, for a weakly coupled cavity (under resonance condition) the main part of 
the field corresponds to a standing wave. Such a field is represented by a sum of 
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real waves, Eq. (8), and a derivation analogous to the outlined above gives 



^'^^tM-^M^I) (16) 

The statistics of eigenstates ip a {v) in a closed system has been rigorously studied 
by Prigodin and co-workers with the help of a zero- dimensional supersymmet- 
ric non-linear q-model|]IO|, [iT| . They studied the joint probability distribution 



P(u\, U2), where Ui = ^(r-j)] 2 (i = 1, 2). For the unitary caseJT0| (broken time- 
reversal symmetry), an expression identical to Eq. (15) (with /j replaced by itj) 



was obtained. For the orthogonal case, Prigodin et al. [11] ended up with a rather 



complicated expression, containing a double integral. Later it has been shown 



by Srednicki [ 12 ] that the expression in Ref.[ll| can be reduced to the function 
given in Eq. (16). He was using the assumption H that a chaotic wavefunction, 
ip a {r), obeys the statistics of a Gaussian random process. This is in complete 
analogy with the standard assumption of the speckle theory |T], [J, as outlined 
above. Again, the difference is that in the speckle theory one usually consid- 
ers propagating waves in an open geometry, whereas Refs. [10-12] study a single 
eigenstate in an isolated system. 

Now, we can analyze the general case of a disordered cavity coupled with arbi- 
trary strength to the external world. The local field is now given by a combination 
of traveling and standing waves, Eq. (10). As a result, the real and imaginary 
parts of the field at two points, ri and r2, are components of a four- dimensional 
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Gaussian vector, <3> T = (Sfa/'ij ^2), with the following covariance ma- 

trix: 



1 



1 + e 2 (1 + e 2 )/ 
e 2 e 2 / 

1 + e 2 )/ 1 + e 2 



e 2 / 
After some lengthy algebra, this leads to 



(17) 



D , T rA exp[-(7 1 + 7 2 )/25] f2*d8 1 d8 2 



25(1 -f) Jo (2tt) 



25(1 - / 2 ) 



x 



(18) 



h cos 201 + I 2 cos 2# 2 + 2/^7^2 ( cos(0x - 2 ) 



costt^ 



Eq. (18) interpolates between a weakly coupled cavity (at the resonance) and an 
open system. 

In conclusion, we consider statistics of radiation in a disordered cavity coupled 
to the environment. The coupling can occur via an opening in the wall of the 
cavity or via absorption in the bulk. For weak coupling, and when the source 
frequency u is close to an eigenfrequency u a , the wave generated in the cavity is 
close to a (standing) eigenmode i/) a ( r ), with only a small admixture of a traveling 
wave. The intensity statistics is defined by the statistics of the eigenfunction 
ip a { T )- F° r strong coupling the system becomesopen and we recover the old 
results of the speckle theory for propagating waves. 

There are many similarities between intensity correlations in open random 
systems and correlations in a single eigenstate of a disordered cavity. There are 
also some differences. It has been already mentioned that, in open systems, the 
field-field correlations decays exponentially for distances p larger than the mean 
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free path I. Therefore correlations described by Eq. (15) are of a short range 
character, p~£. For distances p 3> t a rather different type of correlation, due to 
diffusion, takes over[|l3~ . 



Finally, let us mention that intensity distributions discussed above, such as 
in Eqs. (1), (9) or (13), apply only to the "bulk" of the distributions. Tails of 
the distributions, corresponding to very large or very small values of J, will show 
significant deviations from the above given expressions and will not be universal. 
Indeed, it is well known that, both in open and closed systems, distributions for 
various quantities (conductance, density of states, |iMr)| 2 ) develop log-normal 



tails[14|. This must also be true for the intensity distribution P(I) discussed in 
this paper. For instance, for a point source placed at r = 0, the field ip w (r) is just 
the Green's function G w (0,r) and the intensity is I = \G W (0, r)| 2 . The Green's 
function can be expanded in terms of the eigenf unctions ip Q (r), and the log-normal 
tail of P(\ijj a \ 2 ) are responsible for such tails in the intensity distribution. 
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